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Abstract 

As a typical self-driven many-particle system far from equilibrium, traffic flow exhibits diverse fascinating non-equilibrium 
phenomena, most of which are closely related to traffic flow stability and specifically the growth/dissipation pattern of 
disturbances. However, the traffic theories have been controversial due to a lack of precise traffic data. We have studied 
traffic flow from a new perspective by carrying out large-scale car-following experiment on an open road section, which 
overcomes the intrinsic deficiency of empirical observations. The experiment has shown clearly the nature of car-following, 
which runs against the traditional traffic flow theory. Simulations show that by removing the fundamental notion in the 
traditional car-following models and allowing the traffic state to span a two-dimensional region in velocity-spacing plane, 
the growth pattern of disturbances has changed qualitatively and becomes qualitatively or even quantitatively in consistent 
with that observed in the experiment. 
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Introduction 

Vehicular traffic congestion is a serious problem all over the 
world, which causes huge economic loss and environmental 
pollution. In the last few decades, traffic flow studies have attracted 
wide interests of scientists from various disciplines [1-6]. This is 
because, on the one hand, the studies are practically important to 
mitigate traffic congestion, and on the other hand, traffic flow is a 
self-driven many-particle system far from equilibrium. The 
competitive, nonlinear interactions among vehicles give rise to 
the formation of diverse fascinating phenomena such as boundary 
induced phase transitions, spontaneous formation of jams, 
metastability, hysteresis, phase-separation, and gridlock formation 
[7-1 1]. Therefore, the study of traffic flow offers an opportunity to 
understand various fundamental aspects of nonequilibrium 
systems which are of current interest in statistical physics, such 
as pedestrian traffic and panics [12,13], granular flow [14], 
bacterial colonies [15], swarming and herding of birds and fishes 
[16], traveling waves in an emperor penguin huddle [17], ant trails 
[18], movement of molecular motors along filaments [19], 
dynamic arrest in glasses [20], oscillations of stock market 
dynamics [21], and the "bull-whip" effect of supply chain [22]. 

Most of the nonequilibrium phenomena of uninterrupted traffic 
flow on highways or expressways are closely related to traffic flow 
stability and specifically the growth/ dissipation pattern of distur- 
bances. However, unfortunately, there have been many contro- 
versies in traffic flow research [23-27]. Some traffic engineers 



believe that all traffic jams are induced by bottlenecks [23]. When 
traffic jams occur for no apparent reason, which is known as a 
phantom jam, it is only because nobody has looked far enough to 
find the reason. On the other hand, two decades ago, some 
physicists found that traffic jams can arise completely spontane- 
ously and no bottleneck is necessary [28-30]. Under the right 
conditions small and local fluctuations, which happen all the time 
on the roads, could trigger a non-local traffic congestion that 
persists for hours. Recently, experimental studies on traffic flow on 
a circular road demonstrated that traffic jams can occur on a road 
without bottlenecks [31-33]. 

In the traditional traffic flow theory, a fundamental notion is 
that there is a unique relationship between the flow rate and the 
traffic density (or equivalently between traffic speed and vehicle 
spacing) under steady state conditions [34-37]. A traffic flow 
might be stable, metastable, or unstable. Disturbances in the 
metastable and unstable traffic flows could grow and develop into 
oscillating traffic via a subcritical Hopf bifurcation, like limit cycles 
generated in nonlinear autonomous systems. The traffic flow 
system is thus characterized as a self-excited (autocatalytic) 
oscillator [35]. In contrast, in the three-phase traffic theory 
proposed by Kerner [38-43], it is supposed that the steady state of 
congested traffic occupies a two-dimensional region in the flow 
density plane. Kerner pointed out that there are three traffic 
phases: free flow, synchronized flow, and wide moving jams. 
Usually phase transitions from free flow to synchronized flow 
occur firstly. At the upstream of the synchronized flow, there exists 
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Figure 1. The velocity of the leading car. The driver of the leading car is asked to drive at 50 km/h in Panel A, 40 km/h in Panel B, 30 km/h in 
Panel C, firstly at 15 km/h, then to drive the car without pressing the accelerator pedal (the velocity is about 7 km/h in this case), and then to drive 
the car at 20 km/h in Panel D. Note that for safety reason, the actual velocity of a car is lower than that shown by the speedometer, in particular when 
the velocity is high. The dashed lines are guide for the eyes. 
doi:1 0.1 371 /journal.pone.0094351 .g001 



a pinch region, in which pinch effect induces small jams. The small 
jams propagate towards upstream, grow, merge, and develop into 
wide moving jams. 

The controversy between the three-phase traffic flow theory and 
the traditional one is still on-going, mainly due to a lack of precise 
traffic data [24-27]. Although there are vast amounts of empirical 
data collected by loop detectors, video cameras and floating cars, 
we are unfortunately still not able to know exactly how the traffic 
flow evolves even in very simple scenarios, such as in a platoon of 
passenger cars following each other without overtaking and led by 
a car moving with constant velocity. This is because, on the one 
hand, as claimed by Daganzo et aL, "no empirical studies to date 
describe the complete evolution of a disturbance" and thus we 
cannot "explain the genesis of the disturbances or the rate at 
which they grow" [23]. On the other hand, the empirical 
observations are quite site-specific and contain many confounding 
factors (e.g., geometry, bottleneck strength, traffic flow composi- 
tion), which prevents us from forming a comprehensive under- 
standing of traffic flow evolution. 

Motivated by this problem, we conducted controlled car- 
following experiments to examine the traffic flow theories. 
Different from the previous experiments carried out on a circular 
road [31-33], our experiments concern a platoon of 25 passenger 
cars on a 3.2-km-long open road section. Since the location and 
velocity of each individual car have been recorded by high 



precision GPS devices, the complete evolution of the disturbances 
and their growth rate can be obtained. Moreover, in the 
experiments, the leading car has been asked to move with 
different constant speed so that various traffic flow situations can 
be studied, which goes beyond in-situ observations of traffic flow 
on roads. Therefore, general conclusions about car-following 
behavior are expected. 

Methods 

Ethics Statement 

The study was approved by the institutional review board of 
University of Science and Technology of China. All drivers 
provided written informed consent. 

Experimental setup 

The experiment was carried out on January 19, 2013 on a 
3.2 km stretch of the Chuangxin Avenue in a suburban area in 
Hefei City, China. There are no traffic lights on the road section. 
Since the road is located in a suburban area and has at least three 
lanes in each direction, there is no interference from other vehicles 
that are not part of the experiment. See Videos S 1 and S2 for two 
typical scenarios of the experiment. 

High-precision GPS devices were installed on all of the cars to 
record their locations and velocities every 0.1 second. The 
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Table 1. Sequence of the cars in the two sets of the experiment. 
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measurement errors of the GPS devices were within ± 1 m for 
location and within ± 1 km/h for velocity. Once the experiment 
starts, the driver of the leading car is required to control the 
velocity of the car at certain pre-determined constant values. 
Other drivers in the experiment are required to drive their cars as 
they normally do, but overtaking is not allowed. In two of the runs, 
the driver of the leading car is required to firstly drive the car at 
15 km/h as accurate as he can, then to drive the car without 
pressing the accelerator pedal (the velocity is about 7 km/h in this 
case), and then to drive the car at 20 km/h. In other runs, the 
driver of the leading car is required to accelerate the car to a pre- 
determined velocity and then maintain the car at the constant 
velocity as accurate as he can. The leading car is equipped with a 
cruise control system, which could be switched on when its velocity 
shown by the speedometer reaches 45 km/h. As a result, the 
fluctuation of the leading car is very small when its velocity is set at 
a value greater than or equal to 45 km/h. We note that for safety 
reason, the actual velocity of a car is lower than that shown by the 
speedometer, in particular when the velocity is high. See Figure 1 
for the examples of the velocity of the leading car. When reaching 
the end of the road section, the car platoon decelerates, makes U- 
turn, and stops. When all the cars have stopped, a new run of the 
experiment begins. 

Due to the length limit of the road section, we only analyze the 
experiment data in which the speed of the leading car is below or 
equal to 50 km/h. When the leading car moves with speed larger 



than 50 km/h, the experimental data are too limited after 
removing the data belonging to the start process, the deceleration 
and U-turn process. We have carried out two sets of experiments, 
in which the sequence of the cars has been changed, see Table 1. 

Results 

Car-following behavior 

Figures 2A and 2B show two examples of the evolution of the 
spacing of car No. 1 (which is the second car in the platoon, see 
Table 1) as well as the velocities of the car and its preceding car in 
the first set of experiment. In Figure 2A, the car velocity is nearly 
constant and around 25 km/h, and the velocity difference 
between the car and its preceding one is always small. Before t 
= 174 s, the spacing fluctuates slightly around 15 m. The 
fundamental notion in traditional car-following theory could be 
correct, provided the observations before t = 174 s is the norm. 
However, after t = 174 s, the spacing begins to fluctuate 
significandy and has increased to more than 30 m, even though 
the traveling speed of that car remains nearly constant. Figure 2B 
shows that the spacing also changes significantly between 15 m 
and 40 m when the cars move with velocity around 41 km/h. 
Figures 2C-2F show some more examples of other cars, and 
similar results are observed. 

We would like to point out fluctuations in car-following have 
been discussed previously, see e.g. Refs. [44-46] . For instance, by 
analyzing data from loop detectors as well as data from vehicle 
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Figure 2. Evolution of velocities of a car and its preceding car and their spacing. Panels A and B show Car No. 1 in the 1 st set of 

experiment, Panel C shows Car No. 2 in the 1 st set of experiment, Panel D shows Car No. 23 in the 1 st set of experiment, Panel E shows Car No. 1 1 in 
the 2 nd set of experiment, Panel F shows Car No. 1 2 in the 2 nd set of experiment. See Table 1 for the sequence of the cars in the platoon in the two 
sets of experiment. 

doi:1 0.1 371 /journal.pone.0094351 .g002 



trajectories, Wagner found that the fluctuation in speed difference 
is around ±2 m/s, while the fluctuation in the net time headway 
can be as large as the mean time headway itself [44]. However, 
note that in the empirical data, vehicles change lane frequently. 
Once a vehicle or its preceding one change lane, the vehicle's 
spacing might change significantly. Therefore, as Wagner 
indicates, "unfortunately, that is no guarantee, that the data stem 
in fact from car-following episodes". Therefore, our experimental 
studies have clarified that the car-following process does lead to 
significant fluctuations in spacing. 

This feature clearly contradicts the fundamental assumption 
that there is a unique relationship between vehicle speed and its 
spacing in traditional car-following models. Here are two possible 
explanations of this feature, (i) In certain range of spacing, drivers 
are not so sensitive to the changes in spacing when the velocity 
differences between cars are small. Only when the spacing is large 



(small) enough, will they accelerate (decelerate) to decrease 
(increase) the spacing, (ii) At a given velocity, drivers do not have 
a fixed preferred spacing. Instead they change their preferred 
spacing either intentionally or unintentionally from time to time in 
the driving process, see also Ref. [44] . 

Spatiotemporal features 

Figures 3 A-3E show some typical spatiotemporal patterns of the 
traffic flow. One can observe the stripes structure, which exhibits 
the propagation, growth, dissipation and merge of disturbances. In 
Figures 3A-3D, in which the leading car is not very slow, the 
velocity fluctuation amplitude of the cars in the rear part of the 
platoon has exceeded 20 km/h. In Figure 3D, in which the 
leading car moves with velocity around 1 5 km/h, some cars in the 
rear part of the platoon have occasionally completely stopped. In 
Figure 3E, in which the leading car moves very slowly with 
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Figure 3. Typical spatiotemporal patterns of the platoon traffic. The car speed is shown with different colors (unit km/h) as function of time 
and car number. The left Panels show the experimental results and the right Panels show the simulation results of the 2D ID model. In the simulation, 
the velocity of the leading car is set the same as in the experiment. 
doi:1 0.1 371 /journal.pone.0094351 .g003 



velocity around 7 km/h, the cars in the rear part of the platoon 
will have to completely stop from time to time. 

To better understand the evolution of the disturbances, Figure 4 
presents the standard deviation ff,,of the time series of the velocity 
of each car over all runs of the experiment, in which the leading 
car is asked to move with the same velocity. Note that the data 
belonging to the start process, the deceleration and U-turn process 
have been removed. One can see that when the velocity vi of the 
leading car is low, G v increases almost linearly with car number 
(Figure 4A). With the increase of v h a v becomes increasing in a 
concave way, i.e., the curve of ff,, bends downward and the 
increment rate decreases (Figures 4B-4D). 

Discussion 

Simulations with traditional models 

To study if the observed traffic features of the 25-car platoon 
can be reproduced with traditional car-following theories, we 
choose four traditional car-following models to perform traffic 
simulation. The details of the models and the parameters used are 
as follows. 



Here Ax, = X; — 1 — Xj is the spacing of car i, which follows 
vehicle i — 1. k is sensitivity parameter. £ is noise to represent the 
stochastic factors in traffic flow, which is simply assumed to be 
uniformly distributed random number between — £i and Ci- 

The OV function used is V(Ax) = 1 1.6(tanh(0.086(Ax-25))+ 
0.913) [47]. The unit of length and time are meter and second, 



respectively. The parameters are K '- 
simulation time step At = 0.1s 



1 s _ l, & = 0.2 m/s. The 



The equation of the full velocity difference (FVD) model [48] : 



(hi 
dt 



■- k[V(Axi) - vi] + Uy t - 1 - Vi ) + 1 



Here X is another sensitivity parameter. The OV function and 
the noise £ are the same as in the OV model. The parameters are 
k= 0.32 s~\ ).= 0.4 s~\ £i = 0.2 m/s 2 . The simulation time 
step A? = 0.1s 

• The equation of the intelligent driver (ID) model [49] : 



The equation of the optimal velocity (OV) model [30]: 



dt 



■. K [V(Axt)-v,] + t 



dvi 
dt 
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Figure 4. The standard deviation of the time series of the velocity of each car. The solid lines suggest possible curve fits. 
doi:1 0.1 371 /journal.pone.0094351 .g004 
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Figure 5. Simulation results (solid lines) of the standard deviation. Panel A shows the OV model, Panel B shows the FVD model, Panel C 
shows the ID model, Panel D shows the inertial model. The simulation results are independent of the constant acceleration of the leading car. The 
scatter points are from the experiment. 
doi:1 0.1 371 /journal.pone.0094351 .g005 




0 10 20 30 40 50 
car number 




0 10 20 30 40 50 
car number 




car number car number 



Figure 6. Simulation results (solid lines) of the standard deviation. Panel A shows the 2D OV model, Panel B shows the 2D FVD model, Panel 
C shows the 2D ID model, Panel D shows the 2D inertial model. The scatter points are from the experiment. 
doi:1 0.1 371 /journal.pone.0094351 .g006 
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Figure 7. Simulation results of the spatiotemporal patterns in 
the 2D ID model. The unit of velocities is km/h. Panel A shows V/ = 
65 km/h, Panel B shows v ( = 60 km/h, Panel C shows v t = 55 km/h, 
Panel D shows V; = 50 km/h, Panel E shows V; = 40 km/h, Panel F 
shows v t = 30 km/h, Panel G shows V/ = 15 km/h, Panel H shows V; = 
7 km/h. 

doi:10.1371/journal.pone.0094351.g007 

Here a is maximum acceleration, b is desired deceleration, T is 
desired time headway, s a is the desired gap (bumper-to-bumper 
distance) between two neighboring cars in jam, y max is maximum 
velocity, / is car length, £ is noise as before. The parameters are 
"max = 80 km/h, T = 1.6 s, a = 0.73 m/s 2 , b = 1.67 m/s 2 , ,f 0 
= 2 m, / = 5 m, ft = 0.2 m/s 2 . The simulation time step 
A? = 0.1s 

• The equation of the inertial model [50] : 



dVi ( Vi T + D\ Z 2 (v,— v,-l) 

—r=A \ 1 . J- fcZfVj — v„er) + C 

dt \ Ax,- J 2(Ax,— D) p ' 

Here ^4 is a sensitivity parameter, D is the minimal distance 
between consecutive cars, v per is the permitted velocity, i is a 
constant, Tis the safely time gap constant, £ is noise as before. The 
function £is defined as Z(x) = (|x| + x)/2. The parameters are A 
= 5 m/s 2 , D = 5 m, ^ = 80 km/h, & = 2 s _ 1, T = 2 s, ft = 
0.2 m/s 2 . The simulation time step A? = 0.1s 

Figure 5 shows the curves of the standard deviation of the 
velocities of the cars. In the simulation, the leading car in an 
initially stopped platoon begins to accelerate with a constant 
acceleration until it reaches velocity vi, then it maintains the 
velocity One can see that quantitatively the simulation results 
significantly deviate from the experimental ones. More impor- 
tandy, the simulated curves increase in a convex way, which is 
qualitatively different from the experimental ones. Note that the 



observed differences between the car-following behavior of the 
models and empirical data has also been discussed, see e.g., 
Refs.[44— 46] and references therein. 

We would like to mention that since the velocity of a car is 
bounded by zero speed and the maximum speed, the standard 
deviation ff,,cannot grow infinitely. Thus, the curve of a v cannot 
always grow in a convex way. 

Simulations with new models 

We develop new models by removing the fundamental notion 
that there is a unique relationship between traffic speed and 
vehicle spacing in the steady state in the traditional models. As a 
result, the traffic state of a car could span a two-dimensional (2D) 
region in the velocity spacing plane. Thus, the new models are 
named as 2D models. The details of the 2D models and the 
parameters used are as follows. 

• The 2D OV and FVD models: 

We assume that the relationship between speed and spacing is 
determined by V(Ax) = max(l 1.6(tanh(0.086( m Ax-25))+0.913),0) 
with a parameter m, which changes with time. For simplicity, we 
assume that m is a uniformly distributed random number between 
mi and m 2 , and m changes in the range with rate p (i.e., in each 
time stepA/ in the simulation, m—> m with probability pAt and 
remains unchanged with probability 1- pAt). In the simulations, 
the parameters are p = 0.15 s~ , m\ = 0.8, m 2 — 1.2. Other 
parameters are the same as before. 

• The 2D ID model: 

We suppose that the drivers do not always maintain a constant 
value of desired time headway T. For simplicity, we assume that T 
is a uniformly distributed random number between T x and T 2 , 
and T changes in the range with rate p. In the simulation, the 
parameters are T, = 0.5 s, T 2 = 1.9 s, p = 0.15 s _1 . Other 
parameters are the same as before. 

• The 2D inertial model: 

Similar as in the 2D ID model, we assume that Tis a uniformly 
distributed random number between Ti and T 2 , and Tchanges in 
the range with rate p. The parameters are p = 0.15 s , T\ = 
1.6 s, and T 2 = 2.4 s. Other parameters are the same as before. 

The curves of the standard deviation of the velocities of the cars 
in the 2D models are shown in Figure 6. One can see that 
although the simulation results are still quantitatively different 
from the experimental ones in Figures 6A, 6B, and 6D, the growth 
of the disturbances changes from a convex way into a concave 
way, which is qualitatively consistent with the experimental ones. 
This implies that the feature of traffic flow stability has 
qualitatively changed after removing the fundamental notion. As 
for the 2D ID model, the growth of the disturbances is not only 
qualitatively but also quantitatively in good agreement with the 
experimental ones, as shown in Figure 6C. Figures 3F-3J show the 
spatiotemporal patterns of the velocities simulated from the 2D ID 
model, in which formation of the stripes is very similar to those in 
the experiment. 

Finally we make some predictions with the 2D ID model, which 
need to be validated or invalidated in future work by experimental 
and/ or empirical data. We simulate a large platoon with 600 cars. 
Figure 7 shows the evolution of the velocities of the cars. One can 
see that when V[ is large (e.g., 65 km/h), the traffic flow is quite 
stable. With the decrease of zv, the traffic flow becomes unstable 
and traffic jams appear more and more frequently. Moreover, the 
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locations where traffic jams form become closer and closer to the 
leading car. When V[ is low (e.g., 7 km/h and 15 km/h), high 
speed traffic flow cannot form between neighboring jams. 

Conclusions 

To summarize, we have carried out large scale car-following 
experiments on an open section of a road. The experimental 
evidence clearly runs against the traditional traffic flow theory. 
Simulations with traditional car-following models reproduce 
qualitatively different growth pattern of disturbances from the 
experiment. We have developed new car-following models, in 
which the fundamental notion in traditional traffic flow theory is 
removed and the traffic state of a car could span a two- 
dimensional region in the velocity-spacing plane. The new models 
are able to qualitatively or even quantitatively reproduce the 
growth pattern of disturbances, which implies that the feature of 
traffic flow stability has changed qualitatively. 

We would like to note that due to the limit of road section 
length, the duration of the stationary state is short when the traffic 
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